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ABSTRACT 

We investigate three potential sources of bias in distance estimations made assum- 
ing that a very simple estimator of the baryon acoustic oscillation (BAO) scale pro- 
vides a standard ruler. These are the effects of the non-linear evolution of structure, 
scale-dependent bias and errors in the survey window function estimation. The simple 
estimator used is the peak of the smoothed correlation function, which provides a 
variance in the BAO scale that is close to optimal, if appropriate low-pass filtering 
is applied to the density field. While maximum-likelihood estimators can eliminate 
biases if the form of the systematic error is fully modeled, we estimate the potential 
effects of un- or mi s-modelled systematic errors. Non-linear structure growth using the 
ISmith et' al. (2003) prescription biases the acoustic scale by < 0.3% at z > 1 under the 
correlation-function estimator. The biases due to representative but simplistic models 
of scale-dependent galaxy bias are below 1% at z > 1 for bias behaviour in the realms 
suggested by halo model calculations, which is expected to be below statistical errors 
for a 1000 sq. degs. spectroscopic survey. The distance bias due to a survey window 
function errors is given in a simple closed form and it is shown it has to be kept below 
2% not to bias acoustic scale more than 1% at z = 1, although the actual tolerance 
can be larger depending upon galaxy bias. These biases are comparable to statistical 
errors for ambitious surveys if no correction is made for them. We show that RMS 
photometric zero-point errors (at limiting magnitude 25 mag) below 0.14 mag and 0.01 
mag for redshift z = 1 (red galaxies) and z — 3 (Lyman-break galaxies), respectively, 
are required in order to keep the distance estimator bias below 1%. 

Key words: cosmological parameters - large-scale structure of Universe - dark 
matter - galaxies: statistics 



1 INTRODUCTION 

Acoustic oscillations in the plasma during the pre- 
recombination epoch have b een detected in the cosmic mi- 
crowave background (CMB) (|Spergel et allijlKKi : . Moreover, 
these oscillations should be imprinted on the distribution 
of matter in the Universe and survive until the present 
epoch. The presence of these oscillations was detected as a 
'hump' in the correlation function of Luminous Red Galax- 
ies (LRG) sample in th e Sloan Digital Sky Survey (SDSS) 
|Eisenstein et al]|2005l ). with some evidence in the power 
spectrum of galaxy distributi on in the 2dF Galaxy Red- 
shift Survey (jCole et alj|2005h . These pioneering measure- 
ments reveal the potential of future galaxy surveys, larger 
and deeper than present ones, to measure the sound horizon 
scale as a function of redshift for epochs up to the present 
day. The physics of the plasma era are well understood and 
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let us compute the sound horizon in absolute units (e.g. me- 
ters). BAO measurements may become the best modern in- 
carnation of standard cosmological tests based on a standard 
ruler, for determining the exp ansion history o f the Universe 
and /or in tests for curvature (lBernsteinll2006l) . 

Not much is known, however, about possible system- 
atic effects which could make the measurement of the sound 
horizon more difficult than expected. The purpose of this 
paper is to quantify the bias in the sound horizon scale in- 
troduced by effects like nonlinear evolution, scale-dependent 
bias in tracer objects distribution and observational selec- 
tion (window function) . We compare these levels of potential 
systematic error to statistical uncertainties expected for fu- 
ture hemisphere-scale galaxy surveys. Is the baryon peak in 
the correlation function a robust measure of a distance once 
set in place by physics in the pre-recombination era? 

Quantitative examination of biases in BAO dis- 
tance estimation require that we specify the means by 
which a single distance estimation will be extracted from 
the surveyed galaxy distribution in some redshift range. 
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Figure 1. The nonlinear dark matter power spectrum (upper 
panel) and respective correlation function (lower panel) for our 
fiducial cosmological model at the redshift 2 = 1. Baryon acous- 
tic oscillations are seen as wiggles in the power spectrum and 
the 'hump' at scale 102.1 h _1 Mpc in the correlation function. 
For scales larger than the baryon acoustic peak the correlation 
function becomes negative (dashed line). 



The existing literature on the precision of BAO esti- 
mates considers the distance-es t imation problem in Fourier 
space (Seo & Eisenstein 2003; Blake & Glazcbrook 2003; 
Angulo etalJ 120051 : iBlake et al.1 l200d : ISeo fc Eisenstein! 



20051 ). where the BAO signature is a set of "baryon wig- 
gles" appearing in the galaxy power spectrum at scales 
0.01 h Mpc" 1 < k < 0.5 h Mpc -1 . A model power spectrum, 
uncertain by a scale factor, is fit to the observed galaxy 
power spectrum. In the limit where the galaxy distribution is 
Gaussian and the model power spectrum is exactly specified 
by theory, this yields the maximum-likelihood estimator for 
the characteristic scale and hence an unbiased, minimum- 
variance estimation. If effects such as nonlinear growth or 
galaxy bias are correctly incorporated into the model, then 
the maximum-likelihood techniques still yield optimal unbi- 
ased estimators. It may remain infeasible, however, to model 
all effects exactly or even know their functional forms, so 



we must consider the potential systematic biases from mis- 
modelled physical effects, or those that are not modelled at 
all. Most current BAO analyses or forecasts incorporate ad- 
ditional free parameters, typically in the form of "smooth" 
functions of k that add to or multiply the baseline "wiggly" 
model, in an attempt to mimic any kind of broadband sys - 
tematic effect jSeo fc Eisensteinll2003l ; lTegmark et al.ll2006t ). 
Marginalization over the nuisance smooth-function parame- 
ters yields the single acoustic-scale estimator. 

The Fourier-domain marginalization technique some- 
what obscures the means by which the estimator may end 
up being biased. The situation is much clearer in real 
space, where the linear-era BAO signature is a 3-dimensional 
Green's function comprising a central peak with a spherical 
shell of radius r 3 , the acoustic horizon scale (jEisenstein et al.1 
2006). This clearly produces a single peak at r a in the real- 
space correlation function. We therefore choose to make the 
location of this correlation-function peak our estimator and 
use this very simple estimator to gauge the magnitude of 
biases from unmodelled systematic errors. We will demon- 
strate that this simple estimator is not far from optimal, and 
hence it is useful to quantify its biases. This simple estimator 
incorporates no particular model for non-linearity or galaxy 
biasing; hence by calculating the degree of bias these effects 
generate in this estimator, we have a worst-case estimate. 
Future analysis of the detailed physics will be capable of 
generating corrections that reduce these biases. 

We choose the correlation-function peak to be as model 
independent as possible and avoid the obscuration of fitting 
a smooth function to the power spectrum. The appearance 
of the baryon peak in the correlation function at scales ~ 
100 /i _1 Mpc is difficult to mimic with the small-scale effects 
which can obscure the higher-order wiggles in the power 
spectrum. The power spectrum analysis is appealing due to 
the independence of evolution of Fourier modes (at least in 
the linear regime) , but the correlation- function peak also has 
a simple statistical analysis in the Gaussian limit. 

As this work was completed, lEisenstem" et all (|2006t) 
published an excellent overview of the relation between the 
Fourier- and real-space views of acoustic oscillations in both 
the linear and non-linear regimes. While their work focusses 
on developing an analytical and physical understanding of 
the BAO non-linearities and biasing, we will take the oppo- 
site and less challenging approach of using toy models and 
halo models for non-linear effects, which may not be physi- 
cally accurate, but yield bias estimates that are sufficiently 
quantitative to identify biases that are most in need of fur- 
ther attention. 

The fiducial cosmology is the same as in Blake et al.l 

(2006) to facilitate comparison: flat ACDM with total mat- 
ter density £l m — 0.3 and baryon fraction ft = 15%. The 
present Hubble parameter is taken to be h — 0.7 in units 
of 100 km/s/Mpc, power spectrum normalization as = 1 
and the primordial spectr al index n 3 = 1. We use the linear 
transfer function given by lEisenstein fc Hul (| 1998h . Nonlin- 
ear corrections to the power spect r um a re based on the fit- 
ting formula given bylSmith et"aH (|2003h or the halo model 
l|Coorav fc Shethl 120021 ). The nonlinear fitting formula was 
inferred from results of the pure dark matter simulations 
(the halo model is also usually calibrated with N-body sim- 
ulations) and it is not a priori clear that it is very precise for 
the case of nonzero baryon contribution. The matter power 
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spectrum and respective correlation function for our fidu- 
cial model are shown in Fig. [1] where a 'hump' around the 
scale 100/i~ 1 Mpc is the BAO feature. Changing cosmologi- 
cal parameters make s the 'hump' scale , amplitude and width 
change accordingly (jMatsubaral 120041 ') . A fast and reliable 
transformation between the Fourier and the real space can 
be p erformed by means of the FFTLog routines (|Hamiltonl 
l2OO0h . 

Our analysis is simplified in the sense that we work in 
real space without distinguishing between transverse and ra- 
dial directions and we do not take redshift space distortions 
into account. The former assumption means that both tan- 
gential and radial directions are taken to scale the same with 
the angular-diameter distance Da (z), so information about 
Da(z) might be obtained from either direction. In other 
words we assume that accurate spectroscopic redshifts are 
available for sources with negligible peculiar velocities, and 
that the product of the angular-diameter distance and the 
Hubble parameter Da(z)H(z) is known rather than having 
Da(z) and H(z) measured independently. These assump- 
tions are unrealistic, but should not spoil the rough esti- 
mates of bias from physical effects other than redshift-space 
distortion. 



2 STATISTICAL ERRORS ON THE SOUND 
HORIZON SCALE 

Before addressing the bias in the correlation-function esti- 
mator of r a , we analyse the estimator's statistical errors due 
to the finite survey volume V a (sample variance) and finite 
density of observed objects n g (shot noise). 

Errors associated with sample varian ce scale roughl y 
with the comoving survey volume as V r 1 ^ 2 (|Tegmarkl | 19971) . 
In the "concordance" ACDM model (jSpergel et al.l 12003). 
the comoving volume of space available for observations 
which cover the whole sky grows from 2.4fa _3 Gpc 3 for 
2 < 0.3, through 52/i" 3 Gpc 3 for z < 1 up to 284/i" 3 Gpc 3 
for z < 2.5. For this reason, as well as the degradation of 
the BAO feature by nonlinearities at z < 1, the statistical 
power of B AO surveys degrades s ignificantly at z < 1, al- 
though see lEisenstein et al.l ([2006) for a potential strategy 
to ameliorate the second effect. 

All of the statistical errors on BAO estimators scale 

— 1/2 

with V s . For illustrative purposes we will consider vol- 
umes contained in very wide redshift bins, Az = 0.5. Pro- 
posals exist for near-term surveys to cover ~ 1000 sq. degs. 
of the sky, for which a redshift bin Az = 0.5 centred at z = 1 
gives an observable volume V s = 1.4 ft _3 Gpc 3 , twice as much 
as used by SDSS collaboration to discover the baryon acous- 
tic peak. A highly ambitious goal would be a spectroscopic 
BAO survey covering half of the sky, for a survey volume 
of 28/i _3 Gpc 3 in the z = 1 bin. We will term these the 
"modest-scale" and "hemisphere-scale" survey scenarios. 

The shot noise comes from the finite number of ob- 
jects sampling the distribution of mass. With finite resources 
(telescope time) a survey of the baryon peak has to be bal- 
anced between volume and sampling densi ty to obtain maxi- 
mal si gnal to noise ratio. As e laborated bv lSeo fc Eisensteinl 
(2003) and lMatsubaral l|2004h these considerations typically 
lead to a relation between the optimal number density of 
objects n g and their power spectrum P(k) in the form 



n g P(k) ~ 1. In th e case of BAO we will assume n g P(k = 
0.2/iMpc" 1 ) = 3 jSeo fc Eisensteinl 120031 ) , unless specified 
otherwise. For z = 1 we obtain n g = 2.5 x 10~ 3 /i 3 Mpc -3 
if galaxies are unbiased relative to the mass distribution. In 
sec. 14.11 we will consider galaxy bias and its effect on the 
baryon acoustic peak. 

In order to compute statistical errors in the sound 
horizon scale measurement let us consider the behaviour 
of the galaxy correlation function £ s about its peak at 
r ~ 100/i _1 Mpc. The characteristic scale where the peak 
in the correlation function is observed will be denoted r c 
and is an estimator, possibly biased, of the sound horizon 
scale r a . We allow for smoothing of the observed galaxy dis- 
tribution with a window Wn(k) with a characteristic scale R 
to suppress the small scale power (i.e. mostly the shot noise) 
which contributes to the noise but not to the BAO signal 
in the correlation function. Next, let us compute the deriva- 
tive F(r) of the correlation function as an integral over the 
power spectrum: 



F(r) 



dr 



1 

~2^2 



dkk 3 P(k)WR(k)ji(kr) 



(1) 



where j m {x) is a spherical Bessel function of the m-th order 
The position of the characteristic peak is given by F(r c ) = 
Thus uncertainty in its position is given by 



" \dF/dr\' 

taken at r = r c . From eq. ([1} we obtain the variance of F(r 
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where a%(k) is the variance of a single mo de of the power 
spectrum, given by cr %(k) = (P(k) + l/n g ) 2 (|Feldman et al.l 
1 19941 ; rTegmarkl 1997). Hence, the statistical uncertainty on 
the position of the characteristic scale r c is given by 



(4) 



c [Jdk k*P(k)Wl(k) (j (kr c ) - &h(krc))] 

where a 2 F is given by eq. ((3]). We adopt a Gaussian filter 
of the form Wn(k) — e~ k R ^ 2 , i.e. we presume the galaxy 
density field to be smoothed with a Gaussian of width R 
before measuring the correlation function. 

In Fig. [2] we show statistical errors u(r c )/r c in the 
baryon-peak position, computed by means of eq. @, vs the 
smoothing scale R, for the modest-scale survey at redshift 
z = 1. The statistical errors on the hemisphere survey are 
expected to be 4.5 times smaller. 

The choice of smoothing scale is a trade between beat- 
ing down the shot noise and affecting the BAO signal when 
smoothing is too aggressive — smoothing with R = 7 /i _1 Mpc 
makes the peak disappear into a "knee" in the correlation 
function. A smoothing scale of R = 3 ft Mpc seems to yield 
a signal to noise close to optimal for different shot noise 
contributions (see Fig. [2}, so we adopt this value for the 
remainder of our analysis. 

The smoothing procedure shifts the position of the peak 
in the correlation function, biasing the recovered character- 
istic scale r c . For our assumed R = 3 /i _1 Mpc filter the bias 
is —0.8%. This will not be a concern in the course of this 
work because we are only interested in relative changes of 
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Figure 2. Relative statistical error of the characteristic scale 
Tc in the real space correlation function vs. filtering scale R. 
Note the competing effects of the shot noise subtraction (small R) 
and baryon wiggles removal (large R), both due to the smooth- 
ing. Shown are predictions for the modest fOOO sq. degs. sur- 
vey covering volume of lAh~ 3 Gpc 3 around 2 = 1 (Az = 0.5) 
and spatial galaxy density (solid lines, from top to bottom) 
2.5 x 10~ 3 h 3 Mpc" 3 , 12.5 x 10~ 3 /i 3 Mpc~ 3 , 25 x 1CT 3 h 3 Mpc -3 . 
The dashed line is the sample-variance-limited case. For the full- 
hemisphere survey, all values should scaled with survey volume 



as V, 



-1/2 



lowering errors b y a factor of 4.5. N onlinear evolution 



has been modelled with the ISmith et al. prescription. 



the peak position with redshift. Furthermore, we could in 
practice correct the peak shift using theoretical spectra and 
knowledge of the smoothing procedure. 

The expecte d errors for the aco ustic scale are similar to 
those claimed bv lBlake et"al] |200fJ) . In Table[T]we present a 
concise comparison between expected errors of our the peak- 
of-the-correlation func tion approach (PCF ) and of the power 
spectrum approac h of iBlake et alj (|2006h (B2006). Because 
Bla ke et alj !|2(.IIKi) consider errors in tangential and radial 
acoustic scale separately, we add them in inverse quadrature 
and these combined errors are shown in Table [T] The errors 
are shown for the modest-scale and the hemisphere-scale sur- 
veys around the redshift z = 1 with a depth Az = 0.5. Also, 
two cases of sampling density are considered, one where 
n g = 2.5 x 1CT 3 /i 3 Mpc~ 3 which is implied by the "opti- 
mality" condition n g P(k — 0.2 /i Mpc" 1 ) = 3, the other 
with zero shot noise contribution. We notice that the for- 
mer case equals to almost sample variance limited. The sta- 
tistical errors implied by the PCF method are the same as 
obtained by B2006 in the case of the modest-scale survey 
and only 3% smaller for the hemisphere-scale one, assuming 
n g = 2.5 x 10 -3 /i 3 Mpc -3 . One may infer from this compar- 
ison that the PCF method of the sound horizon estimation 
works well although the estimate of errors is likely too opti- 
mistic (see Sec. [I}' 

Because the peak-of-correlation-function estimator is 
nearly as precise as maximum-likelihood power-spectrum- 
fitting techniques, it is worthwhile to use this simple model- 
independent estimator to investigate bias properties. 



Table 1. Fractional errors on the estimated acoustic scale for 
different surveys obtained be means of the peak-of-the-correlation 
function method ( PCF) described in this paper and the power 
spectrum results of lBlake et al.l (2006) (B2006). The surveys are 
assumed to be around the redshift 2 = 1 and to have the depth 
Az = 0.5. We applied a usual smoothing scale R = 3 h _1 Mpc. 
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Figure 3. The ratio of the nonlinear matter power spectrum 
for a model with baryon fraction = 17%, total matter den- 
sity Q m = 0.27 and ag = 0.9 to the CMB normalised linear 
power spectrum Pjj (fc) with no baryons for redshifts 2 = 3 (long 
dashed), 2 = 1 (short dashed), z = 0.3 (dotted), 2 = (solid). 
The normalised linear power spectrum is also shown (dot dashed) . 
Nonlinear evolution modelling in based on the lSmith et al.l (2003) 
fitting formula and can be compared to the numerical results of 
ISeo fc Eisensteinl d2005h in their Fig. 1. Note that the cosmologi- 
cal model assumed here differs from the one used throughout the 
paper. 



3 EFFECTS OF NONLINEAR EVOLUTION 

The nonlinear evolution of matter perturbations in 
the Universe and its effect on the baryon wiggles in 
the power spectrum have been st udied using numeri- 
cal si m ulations (e.g. Meiksin et al.l (|l999l ): ISpringel et al.l 
(|2005h : ISeo fc Eisensteinl (120051 )1 . Nonlinear evolution cou- 
ples Fourier modes of the matter distribution, suppresses 
power at intermediate scales (k ~ 0.05 h Mpc -1 ), ampli- 
fies power at small scales k > 0.1 h Mpc -1 and partially 
erases the baryon wiggles. This degrades the precision 
of BAO measur e ments for given survey volume, although 
lEisenstein et al.l (|2006h suggest that some of this degrada- 
tion is reversible using corrections for bulk flows. 

Numerical simulations have not, however, been use- 
ful for estimating biases in BAO estimators, because they 
remain too noisy in the range of scales interesting from 
the point of view of baryon wiggles, roughly from k ~ 
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Figure 4. Effect of the nonlinear matter evolution on the mea- 
sured sound horizon scale from the real space correlation function. 
We plot the characteristic scale position r c as a function of red- 
shift for our fiducial cosmological model (solid line) . The horizon- 
tal line (dashed) shows the characteristic scale for the linear mat- 
ter evolution, r c = 102.4 /i -1 Mpc. The error bars represent sta- 
tistical errors in the baryon-peak position from the hemisphere- 
scale survey, assuming that we use galaxies lying in bins of width 
Az = 0.5. 



0.01 /iMptT 1 through k ~ 0.5 ftMpc" 1 jSpringel et all 
120051 ). To estimate the bias in a survey of given volume, 
the simulation volume must be significantly larger than the 
survey volume. Sufficient resources have not yet been avail- 
able to simulate volumes comparable to those to be surveyed 
by a hemisphere-scale survey at 2 > 1, for example. 

One must hence at present rely on analytic arguments 
to investigate the bias of BAO distance estimations. We base 
ou r analysis of nonline ar effects on a fitting formula proposed 
bv lSmith etHI (|2003h which was inspired by the halo model 
and calibrated with N-body numerical simulations. We in- 
clude the effects associated with baryons in the fo rmula via 
the linear transfer function |Eisenstein fc Hulll99cf ). The va- 
lidity of this description of nonlinear evolution has not yet 
been tested extensively, although the formula is able to de- 
scribe reasonably the characteristic features in the power 
spectra of simulations that include baryons. In Fig. [3] we 
show predictions for the power spectrum behaviour given 
by the fitting for mula, which can be com pared to the nu- 
merical results of ISeo fc Eisensteinl (|2005h . Our plot can be 
compared to their Fig. 1; the overall shape of the power 
spectrum is recovered reasonably well, however the wiggles, 
especially those from the third one on, seem to be more 
prominent (less erased) in our description than in the simu- 
lations. 

The influence of nonlinear effects on the measured char- 
acteristic scale r c is shown in Fig. [4] The bias versus the 
linear-regime value of r c — 102.4 h~ 1 Mpc is negative for 
small redshifts and positive for large ones and for sufficiently 
early epochs becomes zero. For epochs earlier than 2 = 1.5, 
the bias in the sound horizon scale is < 0.1% and below the 
statistical error for the hemisphere-scale survey. At redshift 
2 = 1 the fractional bias is Ar c /r c = —0.3%, very close 
to the expected statistical errors from the hemisphere-scale 



survey at that epoch. For redshifts smaller than 2 = 1 the 
nonlinear evolution begins to affect first two baryon wiggles, 
and the bias approaches Ar c /r c — —2.4% at z = 0. In the 
< 2 < 1 range, the bias is comparable to the statistical er- 
rors in the hemisphere-scale survey. Hence the non-linearity 
bias would not greatly dominate the error budget for any 
feasible survey at any redshift range, and in fact would be 
unimportant for any survey that does not approach hemi- 
sphere scale. This is true even without having made any cor- 
rection to the estimator to account for nonlinearities. 

It is worth mentioning that the bias du e to the influence 
of da rk energy on the growth of structure (|McDonald et al.l 
2005) is negligible. 



4 THE GALAXY BIAS AND THE 
BARYON-PEAK POSITION 

The distribution of galaxies is a biased tracer of mass dis- 
tribution and depends on the type of g alaxies considered 
|Zehavi et al.ll2005l ; llvlarinoni et al.ll2005l ). In order to quan- 
tify the effect of galaxy bias on the measured sound hori- 
zon scale, we must specify the scale dependence of the bias 
which, unfortunately, is poorly known at present. One com- 
monly describes the galaxy bias in terms of the multiplica- 
tive factor b(k) relating matter power spectrum P m (k) and 
the galaxy one P g {k) as P g (k) — b 2 (k)P m (k). Obviously any 
bias that is independent of scale leaves the location of the 
peak in the correlation function unchanged. 

4.1 The galaxy bias in the halo model 

The halo model is a statistical description of the structure 
of the Universe based on two basic observations: on large 
scales the matter evolution can be described by perturba- 
tive models and on small s cales matter clusters in to bound 
haloes of a given profile (see lCoorav fc Sh eth (2002) for a re- 
view). The matter and galaxy 2-point functions are divided 
into two-halo and one-halo terms, so the matter and galaxy 
power spectra become 



Pm, g {k) = flin(fc) ( / dMb h (M)n h (M) wZ] g u m , g (k) 



+ 



dMn h (M)w ( ^] g u 2 mig (k) 



(5) 



where is the halo bias, is the halo mass function, 
Um(k) (u g (k)) is the normalised mass (galaxy) distribu- 
tion profile of a halo of mass M in the Fourier space. 
The weight functions for matter and galaxy power spectra 



have the following form for the two-halo term w 

(2) 

Wg 

„.,(!) 



M/p, 

— (N) I n g and for the one- halo term u)„ = M 2 /p 2 , 
= (N(N — 1)) /fig, respectively. N is the number of 
galaxies of a given type occupying a halo of mass M. A 
form of first and second moments of distribution of N must 
be specified and depends on the means of populating haloes 
with galaxies. 

The relation between the mean number of galaxies 
(N) in a given halo and the halo mass M is dubbed the 
halo occupation distribution (HOD). Usually, the HOD for 
luminosity-threshold samples is described by: the minimum 
mass M m i n of a halo in which galaxies may form; the mass 
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Mi of a halo which contain one satellite galaxy on av- 
erage; and the slope a of the relation between number 
of satellite galaxies per halo of a given mass: (N(M)) = 
(1 + (M/M 1 ) a ) H(M - M min ), where H is a step function. 
The details of the halo mode l parameterisation we use can 
be found in iHu fc Jain! l|2004h . In this framework we can ob- 
tain the galaxy bias as a function of scale and galaxy type 
provided that we know the HOD. The specifi c form of an 
HOD was first suggested by theoretical studies |Seliakll2000l; 
Berlind fc WeinbergfeOO^ ; IKravtsov et~aill2004l ; [ Zheng et alj 



20051 ). 

Given the galaxy bias, b(k) = */ Pg{k)/P m {k), or equiv- 
alently the galaxy spectrum P g (fc) , from the halo model, we 
can calculate the correlation function and hence the charac- 
teristic scale r c of its peak. In this halo-model calculation, 
one finds that non-linear growth and biasing does not bias 
the baryon-peak position r c at all (within numerical errors 
which are about 0.1% in r c ) for halo mass thresholds up to 
~ 100M*. This makes sense when one realizes that in the 
halo model the acoustic oscillations are present only in the 
two-halo term, which is assumed to be proportional to the 
linear matter power spectrum, eq. (JSJ). The one- halo term is 
completely featureless at the acoustic scale unless one con- 
siders extremely massive haloes, so the peak of the correla- 
tion function is unchanged. Another way to state this is that 
the halo-model b 2 (k) has small oscillations, coherent with 
the baryon oscillations, that arrange to leave the correlation 
peak unchanged. The scale d ependent galaxy bias fr om the 
halo model was considered bv lSchulz fc Wh ite (2006). Thor- 
ough analysis of the halo bias, based on the perturbation 
theory and numeri cal simulations, was presented recently in 
ISmith et al.l (|2006f h Their approach extends the standard 
halo model by accounting for higher order corrections to the 
linear power spectrum and relaxing a common assumption 
of the linear, scale independent bias in the two-halo term. 
They also discuss potential impact of these improvements 
on the BAOs. 

The near-exact invariance of the correlation func- 
tion peak may be considered a peculi arity of the halo 
model. More sophisticated ha l o models (ISheth et al.l l200ll; 
lYang et al l 120031 ; IZhend 12004 iTinker et all 120051 ) account 
for the effect of previrialisation in the power spectrum which 
the standard model of eq. (0 fails to describe properly. It 
is usually done by replacing the linear matter power spec- 
trum with the nonlinear one in eq. ([5J and accounting for 
halo exclusion in the two-halo term. We do not introduce 
these improvements in our present analysis because the re- 
sultant shifts in the PCF are those we ascribed to non-linear 
growth in Sec. [3] e.g. Fig. 2] In this section we wish to iso- 
late PCF shifts due solely to scale-dependent bias, so we opt 
to describe the nonlinear bias by a generic smooth function, 
(see Sec. I4.2jl that does not share the halo model's excep- 
tional behaviour, and is simply described by the maximum 
and minimum bias values and the fc range over which the 
bias varies. We will use the halo model to determine reason- 
able ranges for these generic parameters, and then estimate 
the shift of the correlation-function peak for the reasonable 
generic functions. 

A three-parameter family of HODs allows one to model 
accurately the projected correlation function of low-redshift 
galaxies (z < 0.07) from the Sloan Digital Sky Survey 
jZehavi et alj|2005h . These authors obtain a good fit to the 



projected correlation function when Mi/M m i n ~ 23, almost 
independent of the luminosity threshold of galaxies in a sam- 
ple. At the same time a rises from 0.9 to 1.2 with luminosity. 
The modelling of the SDSS galaxy correlation function is 
done assuming that the distribution of galaxies within dark 
matter haloes follow the distribution of the matter. 



Unfortunately, such a detailed modelling of an HOD 
is not available for more distant galaxies so we have to 
base our considerations for redshift z>0.1 on results ob- 
tained from num erical simulat i ons. H ydrodynamical simu- 
lations carried by IZheng et al.l (|2005 h broadly support the 
picture which emerges from already me ntioned the SDSS 
galaxy clustering data l|Zehavi et al.ll20"05T ) . The former show 
that both characteristic masses, M m i n and Mi, scale simi- 
larly with galaxy baryonic mass and their ratio for z — 
is Mi/M m i n ~ 14. Although the numeric value of the ratio 
differs from the one obtained from the SDSS galaxy sam- 
ple, the scaling of M m i n and Mi with galaxy baryonic mass 
(simulations) and luminosity (observations) is very similar. 
On the other hand , dark-m atter-only simulations conducted 
by IKravtsov et al.1 (120041 ) show that the relation between 
number of subhaloes and the mass of the host halo is al- 
most linear (a = 1) in a very wide range of halo masses 
(Mi /M m i n > 3) and for cosmic epochs z < 5. Another impor- 
tant result of these simulations is quantifying an evolution of 
the ratio Mi/M m j n with redshift that implies Mi/M m i n ~ 29 
for z = 0, Mi /M mi n ^ 19 for z = 1 and M i /M min ~ 17 for 
z = 3 (see Fig. 5 in IKravtsov et al.l (|2004| y). 



Encouraged by the success of the halo model in describ- 
ing the galaxy clustering we use it for the description of 
the galaxy bias. In this work we consider two HOD models, 
one parameterised by M m i n , Mi and a as mentioned above 
( "central+satellite" scheme) and the other where number of 
galaxies in a given halo is proportional to its mass for haloes 
more mass ive than Mmin (s i milar to the "mass weighted" 
scheme of ISeo fc Eisensteinl ((2005J)). The main differences 
between those two schemes are the following. In the "cen- 
tral+satellite" description if there is at least one galaxy in a 
halo then exactly one galaxy is at the halo centre whereas in 
"mass weighted" scheme all galaxies are distributed within 
a halo. This implies that in the former case the small scale 
power is boosted by the presence of the galaxy in the cen- 
tre, in the latter case small scale power is "smoothed" in a 
way similar to the dark matter distribution. These effects are 
shown in Fig. [5]where we present the galaxy bias at z — 1 as 
a function of scale assuming those two schemes of populat- 
ing haloes with galaxies. The threshold mass M m j n is chosen 
to be 0.1M*, M* or 10M* where M* = 8 X 1O 11 /i" 1 M . We 
note that the bias in the case of "mass weighted" scheme 
shows a rapid change on scales which are potentially rele- 
vant for a behaviour of the baryon wiggles (fc around a few 
tenths of ftMpc -1 ). This is because the most rapid change 
in the bias appears where the galaxy power spectrum is be- 
coming dominated by the one-halo term, whereas the matter 
power spectrum is steep (slope ~ —1.5) in this regime. These 
circumstances appear at different scales, depending on the 
parameterisation of the HOD and the galaxy distribution, 
as seen in Fig. [5] 
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Figure 5. The galaxy bias from the halo model for two schemes 
of populating haloes with galaxies: when number of galaxies oc- 
cupying a halo is proportional to the halo mass (solid lines), and 
when there are central and satellites galaxies with Mi = 20M m i a 
(dotted lines). The minimum mass of the haloes inhabited by 
galaxies for each occupation scheme are 0.1M*, and 10M* 
(from bottom to top for each family of curves). Redshift z = 1 is 
assumed for which M* = 8 X 10 11 1i~ 1 Mq. The bias relation as 
given by the generic function described in Sec. 14.21 is overplotted 
(dashed lines). 



4.2 A toy model 

The halo model confines us to the very specific outcome 
of zero systematic which the Universe may not exhibit. 
We should therefore investigate what possible form of the 
galaxy bias would imply substantial change in the baryon- 
peak position. As mentioned in Sec. 14.11 the scale depen- 
dence of the galaxy bias in the "mass weighted" case as 
shown in Fig. [5] can be described in the interesting range of 
scales (0.01 h Mpc -1 < k < 1 h Mpc -1 ) as a smoothed step- 
like function. A good description can be provided by the 
function of the following form 



logb 2 (fe) = 



log ■ 



erf 



logfc/fc„ 



V2A a 



+ log (6+6-), 



(6) 



where we use the error function and characterise the bias 
b(k) by its asymptotic values at large and small scales, 6_ 
and 6+, respectively (when b+/6_ > 1 we have bias, in the 
other case - antibias). In our analysis it is the ratio b+/b- 
which is important as multiplying the power spectrum by a 
constant factor does not change the position of the peak but 
its amplitude only. The other parameters are the scale k m 
where the most rapid change of the bias occurs, and the loga- 
rithmic interval in k, A CT , centred on k m , over which the bias 
changes by 68.3%. A simple function should suffice if we are 
trying to recover trends in the baryon-peak bias due to the 
galaxy bias and not to obtain its precise values. We exclude 
situations where the bias has a functional form containing 
oscillations of any kind (which was the case in the halo model 
description, see Fig. We expect that the most promi- 
nent effect on the baryon-peak position occurs when the bias 
changes rapidly on scales where baryon wiggles are present 
in the power spectrum, 0.01 foMpc -1 < k < 0.5 ftMpc -1 . 



We now check whether the baryon-peak of the corre- 
lation function is a robust measure of the sound horizon 
in the presence of scale-dependent, smooth deviations in 
the power spectrum. In Fig. [6] we present results of apply- 
ing the galaxy bias of a form given by eq. ([6]) to the lin- 
ear matter power spectrum. The typical scale of the bias 
change is fixed to A CT = 0.4, as suggested by the halo model. 
We take into account only the "mass weighted" case as the 
"central+satellite" scheme introduces negligible bias to the 
acoustic scale. 

The largest effect can be introduced by the galaxy bias 
when the most rapid change in bias takes place around 
the first baryon wiggle in the power spectrum (k m ~ 
0.02 h Mpc -1 ) which itself overlaps the matter power spec- 
trum turnover. The correlation-function peak shift persists 
to k ~ 0.08 h Mpc -1 , around the second peak. It is diffi- 
cult to imagine that the bias has a very rapid change on 
scales as large as the matter power spectrum turnover, which 
is roughly equal to the largest scale of possible physical 
interactions in the pre-recombination Universe. For scales 
km ~ 0.5 h Mpc -1 , where it is more reasonable to consider 
scale dependence of the galaxy bias, a very strong bias ampli- 
tude of order several is required to change the baryon-peak 
position by more than 2%. If we look to the halo model 
for an estimate of the amplitude of scale-dependent bias, 
we find the shift in the characteristic scale r c is below 1% 
at 2 = 1 for any threshold mass below 3M* (the survey 
depth is Az = 0.5). The halo-model parameters for our 
toy model further suggest that, for surveys of galaxies in- 
habitating haloes of mass M* a bias of the acoustic peak 
is 0.7% at the redshift z — 1, whereas expected statisti- 
cal error for the modest-scale survey is 1.5%. For redshifts 
z — 2 and z — 3 respective systematic errors are 0.5% and 
0.2% and statistical ones 1% and 0.95%. Statistical errors for 
the hemisphere-scale survey are smaller by a factor of 4.5. 
The bias becomes comparable to statistical errors for the 
hemisphere-scale survey, but not for the modest-scale one. 
This is without any correction to the estimator r c . Numeri- 
cal models should eventually allow for corrections to reduce 
the bias below statistical errors ever for hemisphere-scale 
surveys. 

It is difficult to assign the minimum halo mass that will 
correspond to future BAO galaxy surveys at z>l. It will 
depend on the selection criteria for galaxies in the future 
surveys. Haloes of mass > M t ar e possible hosts of SDSS 
Luminous Red Galaxies at z ~ 0.3 (|Zehavi et al . 2005). Sim- 
ilar studies of halo occupancy have yet to be done at higher 
z, but the appropriate halo-model mass threshold M m i n will 
likely be in the range 0.1 — 10Af*. 

When estimating statistical errors we assumed no 
galaxy bias. In practice, the bias is likely to help in low- 
ering statistical errors either through allowing for a denser 
sampling or covering larger area while using the same re- 
sources. 



5 SURVEY WINDOW FUNCTION 

Another effect which requires attention is the possible in- 
fluence of the survey selection function on the baryon-peak 
position in the correlation function. Let us assume that the 
observers of our fictitious real-space survey have estimated 
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Figure 6. Contour plot of the fractional bias of the baryon- 
peak position as a function of scale k m at which the most rapid 
change of the galaxy bias appears and the relative amplitude of 

the galaxy bias 6_|_ /b . The range of scales where bias changes 

the most rapidly is fixed to A CT = 0.4. Open squares mark values 
of parameters which describe by means of the smooth function 
the galaxy bias obtained from the halo model (see Fig. [5j assum- 
ing minimal halo masses of 0.1M*, M*, 10M* (from bottom to 
top). Shaded regions show parts of the parameter space where 
statistical errors on baryon-peak position dominate over the sys- 
tematic shift. Statistical errors are assumed to be for two surveys: 
the modest-scale (light shaded region) and the hemisphere-scale 
(dark shaded region). 



a spatial selection function Wo(r), but the correct selection 
function differs by a small multiplicative correction Sw(r). 
The observed galaxy density will be 



Pg (r) = n g W (r) (1 + S w {r)) (1 + i true (f)) . 



(7) 



The experimenters will estimate the overdensity using their 
estimated selection function, so will calculate an (incorrect) 
local overdensity 



(1 + S ohs (f)) = (1 + 8 w {r)) (1 + S tIue (r)) . 



(8) 



The observers' calculated correlation function £ b s follows: 
1 + U s (r) = (1 + VM) (1 + £tru.(r)) , (9) 

where we assume that the selection function is uncorrelated 
with the true density field, and £w is the autocorrelation of 
the fractional errors &w in the selection-function estimate. 

What level of residuals in the window function are al- 
lowed in order not to bias the sound horizon scale determina- 
tion more than predicted statistical errors? We are interested 
in scales which are close to the baryon peak in the correla- 
tion function, so we can approximate correlation function 
£true(r) to second order about the true acoustic scale 



ftrue(Y) = £ truo (r s ) ( 1 - 



(r- 



1ol 



(10) 



ctq has an interpretation as the width of the baryon peak. 
The characteristic scale r c is measured from the position of 
the peak in the observed function £ b s (r), so the derivative 
of the £ bs is zero at r c . Hence the shift in the baryon-peak 



position due to the presence of spatially-varying error in the 
selection function is 



pn(l + Mr c ))]' 



(1 + £truc(r s )) 
Ct'rue ( r s ) 



(11) 



where the prime means the derivative with respect to r and 
only linear terms in Sr 3 are retained in the expression for 
Ctruc- Note the bias depends only upon the properties of £w 
at the acoustic scale. 

From eqs. (|10[1 and follows that for the power law 
£w( r ) with an exponent —a and an amplitude £w *C 1 in 
the vicinity of the observed peak we obtain to the leading 
order 



5r a ( go 

— = -ai- 
rs \ r s 



jw{r c ) 

£true (^s ) 



(12) 



The narrower the acoustic peak, the less easily it is moved; 
and the shallower (r) is at the acoustic scale, the smaller 
a bias in the peak position is produced. 

We can make an estimate of the level of systematic er- 
rors due to the unknown properties of the window func- 
tion for a given survey. As an example let us consider 
the survey at the redshift z ~ 1 in which case we had 
r s = 102.1 /i _1 Mpc, £true(r s ) = 6.33 x 10" 4 and a = 
10.2/i -1 Mpc where we assume no galaxy bias (Sec. [2j and 
a = 2. Thus eq. (|12|l implies that to be able to beat down 
the systematic error to the level of 8r s /r s , e.g. 1%, the cor- 
relation function of the window function correction cannot 
be larger than 



\Zw(r c )\ 



3.2 x 10" 



(?r 



£true (f"s ) 

6.33 x 10- 



102.1 fe-'Mpc 



lO^ft-iMpc 



Sr s /r s 



1% 



(13) 



Hence, eq. (|13[1 shows that unknown variation in the selec- 
tion function on scales of the baryon peak should be < 1.8% 
if we want to bias our results no more than 1%. For larger 
redshifts the allowed variation would have to be smaller — 
1.3% and 1% for z — 2 and 3, respectively. Scaling of the 
variation with redshift is mainly due to the linear suppres- 
sion of the mass fluctuations and narrowing of the baryon 
peak (several times smaller effect). A survey of biased galax- 
ies makes £true(r s ) rise, allowing for a larger error in the 
selection function without biasing r c by more than the sta- 
tistical error. 

An example of an effect which can bring about correla- 
tions on large scales considered above is the extinction pat- 
tern on the sky wh ich is known to be c orrelated on scales 
of tens of degrees (jSchlegel et al.|[l998l ). Also, for surveys 
around redshift z — 1(3) the acoustic scale is about 2.5(1.3) 
deg (as opposed to ~ 6 deg for the SDSS LRG sample). This 
is around (or less) than a size of a field of view of planned 
wide-field surveys, thus a different magnitude limit for each 
telescope po inting may introduce cor relations on the inter- 
esting scale (|Guzik fc Bernsteinll2005l ). Moreover, significant 
uncertainty to a measured window function may be intro- 
duced by a complicate d geometry of a survey |Miller et al.l 
120021 : ICole et alJ[2005h . 

We investigate a model in which each telescope point- 
ing has an independent magnitude calibration error, and 
relate the RMS magnitude error to the systematic bias on 
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the baryon peak position through the survey window func- 
tion as assumed in eq. (JHJ. The correlation function of the 
angular window functio n for calibration errors i n circular 
fields of diameter 6 a is (|Guzik fc Bernstein! 1 2005) £ w {6) = 
2/ttE 2 (arccos (<9/(9 a ) - (9/8 a ) - (0/0 tt ) 2 ) if 9 < 6 a and 
zero otherwise. Now let us compute the variance, E 2 , which 
is given as a fractional uncertainty in the density of objects 
in the real space due to the survey limiting magnitude error 



4>(L/ Li 



. dh 
' dm 



5(n g W ) = 

n g W f™ Lt c?{L/L*)dL/L+ 



5m, 



(14) 



where (j)(L) is the luminosity function in the Schechter form 
4>{L)dL — (f>*(L/L*) a exp(— L/L*)d(i/L*) and mi im is a 
survey limiting magnitude. For a s urvey of red (blue) galax- 
ies at z ~ 1.1 jFaber et al.l 12005) the luminosity function 
is described by a characteristic density of = 1.5(8.4) x 
10 _3 ( /i -1 Mpc) -3 , characteristic absolute magnitude M* = 
-21.58(-21.25) and the faint-end slope a = - 0.5(-1.3). At 
z ~ 3, observations of Lyman-break galaxies (|Steidel et alj 



1 19991 ) yield 0* = 48.0 x 10~ 3 ( /i^Mpc)" 3 , M* = -21.68, 
a = —1.6. Also, let us assume that the survey telescope 
aperture is 9 a = 3 deg and the limiting apparent magni- 
tude is 25 mag. A predicted position of the baryon peak 
at z = 1(3) (see Sec. [TJ is 102.1(102.5) /i _1 Mpc, its width 
10.2(9.5) fe _1 Mpc and the matter correlation function ampli- 
tude at the peak scale 6.33(3.16) x 10 -4 . Then from eq. (|ll[) 
we obtain that the magnitude calibration error 5m should 
be smaller than 0.14(0.05) x b 2 for survey of red (blue) 
galaxies at z — 1 if we require an error on the baryon 
peak position be smaller than 1%. Note that the calibra- 
tion depends on the galaxy bias b. Respective calibration 
er ror for z ~ 3 yields 0.01 x b 2 . If we extrapolate results 
of iPadmanabhan et al.1 12006) to the bias of Luminous Red 
Galaxies to z — 1 we could expect the bias of this type of 
galaxie s is b ~ 2. Blue gala xies are expected to be much less 
biased l|Zehavi et alj 2005). Thus, the magnitude calibration 
for the redshift survey of the red galaxy sample at z ~ 1 
should not present any practical challenge. The considered 
effect can have more impact on the blue-galaxy sample due 
to the steeper faint end of their luminosity function. Also, 
for a survey at z ~ 3, the required calibration accuracy may 
be more diffic ult to obtain, although a large bias, of order a 
few o r more (|Porciani fc Giavaliscd |2002| ; [Kashikaw a et all 
2006), can help. Predictions for z ~ 3 can vary because the 
limiting magnitude is on the exponential part of the lumi- 
nosity function. 



6 CONCLUSIONS 

We have considered three types of possible systematic effects 
in the estimation of the cosmological distances by means of 
the baryon acoustic peak in the correlation function. Our 
analysis accounts for the nonlinear evolution of structure, 
the galaxy biasing with respect to the underlying mass and 
the effect of imperfect recovery of the survey selection func- 
tion. In each of these cases we compare systematic to statis- 
tical errors that one could expect for a 1000 sq. degs. and a 
half of the sky spectroscopic surveys. 

The nonlinear evolution biases the estimated acoustic 
scale for redshift z < 1.5. For redshift z — 1 (and the depth 
of the survey Az = 0.5) the systematic error is 0.3% and 



grows up to 2.4% at z = 0, close to the predicted statistical 
errors for the hemisphere-scale survey. Statistical errors are 
expected to dominate the systematic ones for any survey 
smaller than this, e.g. the modest-scale survey yields statis- 
tical errors larger by a factor of 4.5 than the hemisphere- 
scale one. Thus, the effect of the nonlinear evolution on the 
measured acoustic scale is expected to be below statisti- 
cal errors for z > 1 in case of planned surveys even with- 
out applying correction for nonline ar evolution. Correctio n 
techniques, like the one proposed bv lEisenstein et all (2006), 
may be useful in unbiasing the acoustic scale even for shal- 
low surveys, 2 < 1. 

The measured acoustic scale may also be affected by 
the scale-dependent galaxy bias which is poorly known at 
present. We quantified a possible effect of the galaxy bias 
based on the toy model supported by the halo model results. 
We showed that for the case of the galaxy bias predicted by 
the halo model it is unlikely to introduce systematic effect on 
acoustic scale exceeding 1% when haloes of mass < 3M* at 
the redshift z = 1 are considered. For less massive haloes the 
acoustic scale bias is smaller. In fact, the bias for M* haloes 
at z = 1 is 0.7% which is lower by a factor of 2 than sta- 
tistical errors for the modest-scale survey (1.5%) and higher 
than those expected for the hemisphere-scale survey (0.3%). 
For higher redshifts, the acoustic scale bias for haloes 
is expected to be at least a factor of 2.5(4.75) smaller than 
the statistical errors for the redshift z = 2(3) in case of the 
modest-scale survey. On the other hand, the systematic er- 
ror could be larger if the change of the bias with scale was 
more rapid or took place at larger scales than suggested by 
the halo model. 

Also, we considered the bias in the acoustic scale mea- 
surement due to the correlated errors in the survey win- 
dow function. In our analysis the systematic error depends 
linearly upon the local slope of the correlation function 
of the selection function residuals and varying inversely 
with the local curvature of the true correlation function. 
Our analysis implies that in order to achieve a bias < 1% 
one has to measure the window function (more precisely, 
£w(~ 100 /i _1 Mpc)) with accuracy better than 2% at z = 1, 
and more accurately for higher redshifts. The constraint can 
be relaxed if the galaxy bias on scales of the acoustic peak is 
greater than 1. Specific results depend on details of a survey 
and a data reduction methods. We considered one model, 
namely independent photometric zero-point errors at each 
telescope pointing. RMS zero-point errors (at limiting mag- 
nitude 25 mag) below 0.14 mag and 0.01 mag for redshift 
z — 1 (red galaxies) and z = 3 (Lyman-break galaxies), 
respectively, are required in order to limit systematic er- 
ror in the peak-of-the-correlation-function estimator to less 
than 1%. In order to gain a confidence in unbiasedness of 
the acoustic scale measurement, techniques to diagnose win- 
dow function problem s have to be developed and applied 
l|Tegmark et al.ll2002l ). 

We compared the variance of the peak-of-the- 
correla tion function esti mator to the power spectrum esti- 
mator (|Blake et al.ll2006l ) and found the former yields virtu- 
ally the same statistical inaccuracies as the latter. This sug- 
gests that our model-independent estimator might be useful 
in measuring of the sound horizon scale in spite of the sim- 
plifying assumptions we made in the course of its analysis. 
Even though its variance may turn out to be larger, its sim- 
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plicity and model independence are sufficient reasons for the 
peak-of-the-correlation method to be practically important. 

In this paper we considered only the real space correla- 
tion function. In the redshift space, a peak in the correlation 
function becomes an elliptical ridge and one has to properly 
modify the peak-of-the-correlation function estimator to be- 
come applicable in this case. Fortuna tely, the large scale red- 
shift space distortions (|Kaiserlll987h do not lead to a bias 
of the baryon-peak position because the relation between 
an azimuthally averaged redshift space power spectrum and 
the real space one is scale independent. How ever, an im- 
proved analysis of redshift space distortions l|Seliakll200i 
IScoccimarrol I2004T ) may lead to different implications. We 
defer an analysis of the redshift space effects to the future 
work. 



ACKNOWLEDGMENTS 

We would like to thank Bhuvnesh Jain, Mike Jarvis, Hee- 
Jong Seo and Ravi Sheth for numerous helpful discussions. 
This work is supported by NSF grant AST-0236702, Depart- 
ment of Energy grant DOE-DE-FG02-95ER40893, NASA 
grant BEFS-04-0014-0018 and Polish State Committee for 
Scientific Research grant 1P03D01226. 



REFERENCES 

Angulo R., Baugh C. M., Frenk C. S., Bower R. G., Jenkins 

A., Morris S. L., 2005, MNRAS, 362, L25 
Berlind A. A., Weinberg D. H., 2002, ApJ, 575, 587 
Bernstein G., 2006, ApJ, 637, 598 
Blake C, Glazebrook K., 2003, ApJ, 594, 665 
Blake C, Parkinson D., Bassett B., Glazebrook K., Kunz 

M., Nichol R. C, 2006, MNRAS, 365, 255 
Cole S. et al., 2005, MNRAS, 362, 505 
Cooray A., Sheth R., 2002, Phys. Rep., 372, 1 
Eisenstein D. J., Hu W., 1998, ApJ, 496, 605 
Eisenstein D. J ., Seo H.-j., White M., 2006, 

|astro-ph/0604361| 

Eisenstein D. J. et al., 2 005, ApJ, 633, 560 
Faber S. M. et al., 2005, |asta>ph/0506044 
Feldman H. A., Kaiser N., Peacock J. A., 1994, ApJ, 426, 
23 

Guzik J., Bernstein G., 2005, Phys. Rev. D, 72, 043503 

Hamilton A. J. S., 2000, MNRAS, 312, 257 

Hu W., Jain B., 2004, Phys. Rev. D, 70, 043009 

Kaiser N., 1987, MNRAS, 227, 1 

Kashikawa N. et al., 2006, ApJ, 637, 631 

Kravtsov A. V., Berlind A. A., Wechsler R. H., Klypin 

A. A., Gottlober S., Allgood B., Primack J. R., 2004, ApJ, 

609, 35 

Marinoni C. et al., 2005, A&A, 442, 801 

Matsubara T., 2004, ApJ, 615, 573 

McDonald P., Trac H., Contaldi C, 2005 , |astro-ph /0505565 1 
Meiksin A., White M., Peacock J. A., 1999, MNRAS, 304, 
851 

Miller C. J., Nichol R. C, Che n X., 2002, ApJ, 5 79, 483 
Padmanabhan N. et al., 2006, | astro-ph/0605302] 
Porciani C, Giavalisco M., 2002, ApJ, 565, 24 



Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Schulz A. E., White M., 2006, Astroparticle Physics, 25, 
172 

Scoccimarro R., 2004, Phys. Rev. D, 70, 083007 

Seljak U., 2000, MNRAS, 318, 203 

Seljak U., 2001, MNRAS, 325, 1359 

Seo H.-J., Eisenstein D. J., 2003, ApJ, 598, 720 

Seo H.-J., Eisenstein D. J., 2005, ApJ, 633, 575 

Sheth R. K., Hui L., Diaferio A., Scoccimarro R., 2001, 

MNRAS, 325, 1288 
Smith R. E. et al., 2003, MNRAS, 341, 1311 
Smith R. E., Sc occimarro R., Sheth R. K., 2006, 

|astro-ph/0609547| 

Spergel D. N. et al., 2003, ApJS, 148, 175 
Springel V. et al., 2005, Nat, 435, 629 

Steidel C. C, Adelberger K. L., Giavalisco M., Dickinson 

M., Pettini M., 1999, ApJ, 519, 1 
Tegmark M., 1997, Physi cal Review Letters , 79, 3806 
Tegmark M. et al., 2006, j astro-ph/0608632| 
Tegmark M., Hamilton A. J. S., Xu Y., 2002, MNRAS, 335, 

887 

Tinker J. L., Weinberg D. H., Zheng Z., Zehavi I., 2005, 
ApJ, 631, 41 

Yang X., Mo H. J., van den Bosch F. C, 2003, MNRAS, 

339, 1057 
Zehavi I. et al., 2005, ApJ, 630, 1 
Zheng Z., 2004, ApJ, 610, 61 
Zheng Z. et al., 2005, ApJ, 633, 791 



© 0000 RAS, MNRAS 000, 000-000 



